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Radiative Upsilon Decay at the Endpoint 
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The standard NRQCD power counting breaks down and the OPE gives rise to color-octet shape functions at 
the upper endpoint of the photon energy spectrum in radiative T decay. Also in this kinematic regime, large 
Sudakov logarithms appear in the octet Wilson coefficients, ruining the perturbative expansion. Using SCET, 
the octet shape functions arise naturally and the Sudakov logarithms can be summed using the renormalization 
group equations. We derive an expression for the resummed octet energy spectrum. 
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The decay of the T in the endpoint region is 
interesting for a number of reasons. Firstly, it is 
a way to test Non-Relativistic QCD (NRQCD) 
in decays, which could be important for extrac- 
tions of a s 0. Secondly, in the endpoint region, 
it is necessary to introduce another effective field 
theory, Soft-Collinear Effective Theory (SCET) 
IU, due to the jet of collinear particle that is pro- 
duced in this kinematic region. Finally, CLEO 
will obtain much more precise data on the first 
three T resonances, so it is timely to see if we 
can predict the decay spectrum. 

The T(15) decays radiatively about 3% of the 
time. When the photon is emitted with maximum 
energy, z = 2E^/My — * 1, it is back-to-back with 
a jet of particles. It is this kinematic region that 
we are interested in predicting Before dis- 
cussing this region in particular, it will be useful 
to review how the calculation is done in general. 

To calculate the inclusive photon spectrum, it 
is advantageous to use NRQCD. The effective 
field theory (EFT) is based on the fact that the 
b quark is heavy, m,b ^> Aqcd- This implies that 
the relative velocity of the the heavy quarks in- 
side the bound state is small, v <C 1. For T, 
v 2 w 0.1. It is therefore beneficial to do a double 
expansion in v and a s . Furthermore, it is possi- 
ble to use NRQCD to prove the short and long 
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distance physics factorize in T decay |J . 
The decay rate in NRQCD is written as 



dT 

dz 



J2C n (z)(0(n)). 



(1) 



The C n are Wilson coefficients, calculable in 
perturbation theory, while the (0(n)) are non- 
perturbative NRQCD matrix elements (MEs). 
The MEs schematically are 



(O(n)) = <T|^ t r"Wr 7 V|T}, 



(2) 



where the T n can contain derivatives, color and 
spin matrices, which can be classified in spectro- 
scopic notation and as a color-singlet or color- 
octet. For instance (O 8 ( 1 5 )) is a ME where T n 
contains only the color-matrix T a . 

The series in Eq. ([!]) is infinite, so to have any 
predictive power, we need to truncate. This is 
possible using the velocity scalings of the MEs. 
Each ME scales as a certain power of v, depend- 
ing on what terms from the NRQCD Lagrangian 
needs to be inserted to have a non-zero overlap. 
For the endpoint spectrum, the relevant MEs are 



(Oi( 3 5x)) 
(OsCSo)) 



v 4 . 



(3) 
(4) 
(5) 



The color-singlet ME (Oi( 3 S*i)) can be related to 
the wavefunction at the origin, and the decay rate 
through this channel is the result obtained in the 
Color Singlet Model f§. 
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At first glance, it appears that the color-octet 
contributions to the rate are tiny, v A « 0.01, 
compared to the color-singlet. However, that is 
not true. To see why, we need to compare the 
rate for each channel, including the Wilson coef- 
ficients. For the color-singlet rate, the T decays 
to a photon and two gluons, so c[ a \ 3 S±) cx a 2 . 
For the color-octet channels, on the other hand, 
the final decay products are a photon and one 
gluon. So Cf ] ^S ) ~ C 8 ( 3 F ) oc ira s , where the 
7r comes from there being one less particle in the 
final state. The color-octet is enhanced perturba- 
tively by a factor of ira s — O(10). 

But that is not all. Since there are only 
two particles in the final state for the color- 
octet decay, the rate is peaked at the endpoint, 
C^CSo) ~ C^ 0) ( 3 P ) cx 6(1 - z). If we com- 
pare the integrated rate in the endpoint region, 
(1 - v 2 < z < 1), we have r5 nd ( 3 S"i) cx a 2 s v 2 for 
the color-singlet, and r™ d ( 1 S' ) - r™ d ( 1 5 ) cx 
ira s v 4 ss oc 2 v 2 for the octet, using the fact that 
numerically v 2 w a s /ir. So in the endpoint region 
the singlet and octet rates contribute equally. 

If we went to higher order in the velocity ex- 
pansion, we would find further contributions, for- 
mally suppressed by higher powers of v, that con- 
tribute equally in the endpoint region. This is due 
to a breakdown of the non-perturbative expan- 
sion. The solution has been known for some time 
Q. The series needs to be reordered as a twist 
expansion, similar to what is done in B decays 
M . The rate is then written as 



pear in the Wilson coefficients at NLO of the form 



dT 

dz 



dk + C n (z,k + )f n (k + )(0(n)), 



(6) 



where f n (k + ) are shape functions, which measure 
the probability for the quark pair (in state n) to 
have light-cone momentum k + . As these are non- 
perturbative functions, they must be modeled. 

We also must worry about the perturbative se- 
ries. The fact that the color-octet rate begins as 
a 6 function is already worrisome. Higher order 
corrections could lead to further problems. This 
in fact is true. The next-to-leading order (NLO) 
perturbative corrections have been calculated for 
the color-singlet (numerically) Q and color-octet 
||. The NLO color-octet rate is singular at the 
endpoint. In particular, Sudakov logarithms ap- 
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(7) 



which become large as z — > 1. If we again inte- 
grate over the endpoint, (1 — v 2 < z < 1), these 
terms give contributions to the rate a 2 log 2 (v 2 ) 
and a 2 log(t> 2 ). Both are 0(a e ), since a s log 2 (v 2 ) 
and a s log(w 2 ) are 0(1). This is an indication 
that the perturbative series breaks down in the 
endpoint. If wc looked at higher order perturba- 
tive corrections, we would get terms of the form 
a s [o!" log 2n (v 2 )], which are also 0(a s ). 

The break down of the perturbative and non- 
perturbative series are due to the same problem: 
NRQCD does not contain all of the correct de- 
grees of freedom for the endpoint region. It is 
missing collinear modes. To correctly describe 
the physics at the endpoint we need to couple 
NRQCD to an EFT that contains these missing 
fields, SCET @. 

The invariant mass of the hadronic jet in the 
endpoint region, M| = (1 — z)M T = A 2 M T , is 
much larger than the energy in the jet, (Ex = 
zM-r/2). We therefore have a multiscale system, 
and EFT techniques are useful to separate the 
scales. In this case, the expansion parameter A = 
VI — z. There are three classes of particles that 
are included in the EFT, depending on the scaling 
of the momentum. They are massless quarks and 
gluons with: 



Collinear 

Soft 
Ultrasoft 



E + p A ~ My, E-p A 
p 1 ' 2 ~ AM T , 
p v ~ AM T , 



A 2 My 



^ - A 2 My. 



The scales in the problem are the hard scale, 
My, the collinear scale, AMy, and the soft scale, 
A 2 A/y- The plan of attack is to integrate out the 
hard scale by matching to the EFT and then run 
to the collinear scale using the renormalization 
group equations (RGEs) of SCET. At this point 
all collinear particles are far off shell and should 
be integrated out. This is done by matching onto 
the soft theory, which is same as the large energy 
effective theory JlO|]. We will then run from the 
collinear scale to the soft scale, at which point all 
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Figure 1. Matching onto SCET. The diagram on 
the left is calculated in full QCD, and then ex- 
panded in A and v. This is matched onto an op- 
erator in the EFT. 



the large logarithms will be in the coefficient func- 
tions. The logs in this case are log A — > log(l — z), 
so by using the EFT RGEs, we have resummed 
all the logs of (1 — z), i.e. the Sudakov logs j|. 

To follow this program, we first need to match 
onto SCET. This is done by calculating the 
graphs in QCD and expanding in A (and in v 
which will also match onto NRQCD). An exam- 
ple of this is shown in Fig. [|. The matching, 
at leading order in A, gives two operators in the 
EFT, for the octet 1 S and 3 P channels. We 
do not get any other operators at leading order, 
which means there are no leading Sudakov logs 
in the color-singlet channel Eu . There are other 
operators suppressed by A |12[ . 

We next need to run using the SCET RGEs. 
This requires the anomalous dimensions of the 
operators. We therefore need to calculate the one- 
loop diagrams in Fig. ^. The anomalous dimen- 
sions of the 1 So and 3 Po operators are the same, 
which is related to the fact that the Sudakov logs 
are the same for both channels |J. Please sec [|J 
for the details. 

We next match onto the soft theory, by remov- 
ing all collinear particles. However, we have a 
collinear gluon in the final state, so to remove 
it, we perform an operator product expansion 
(OPE). The result is a non-local operator, sep- 
arated along the light-cone. This is shown dia- 
grammatically in Fig. [| The operator we match 



Figure 2. The one- loop diagrams needed to calcu- 
late the anomalous dimensions of the SCET op- 
erators. 



onto is of the form 

O n (x) = ^T /n x S(l -x + iD+/My)x^r n ^, (8) 

where the derivative is in the + light-like direc- 
tion. The rate is now 



dT 

dz 



J dx^2c n (x - z)f n (x){0(n)), 



(9) 



where f n (x) — {O n {x)) / (0(n)) . This is just 
the rate including the shape function. By using 
SCET, the shape function appears naturally. 

We now run down to the soft scale. For the 
details, see H. At this point, there are no large 
logs in the operators. We have summed all the 
Sudakov logs into the Wilson coefficients. To plot 
the result, we use the model of Jr|], introduced for 
B decays. For T, we need to give the shape func- 
tion some unknown first moment of order Aqcd 
M. The results are shown in Fig. ^. The dashed 
curve is the resummation without the shape func- 
tion. The dotted curve is the singular terms in the 






Figure 3. Matching onto the soft theory, by doing 
an OPE. 
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Figure 4. The differential decay spectra near the 
endpoint region in arbitrary units. The curves 
are described in the text. 



Figure 5. The different channels for the spec- 
tra near the endpoint region. The curves are de- 
scribed in the text. 



one-loop result, and the dot-dashed curve is these 
terms convoluted with the shape function. The 
solid curve is the resummation convoluted with 
the shape function. Note that the shape func- 
tion without resummation and the resummation 
without the shape function give similar results. 
However, both are necessary. 

At this point we compare to the color-singlet 
results. To do this we need the color-octet MEs. 
In Fig. the solid curve is the color-singlet result. 
The dotted curve is the resummation without the 
shape function, the dashed is resummation with 
the shape function. The two sets of curves corre- 
spond to scaling the octet MEs from the singlet 
by factors of t> 4 /10 and u 4 /100. As can be seen, 
to get comparable results, the naive scaling (v 4 ) 
seems to be off by a factor of 100 

However, before a meaningful comparison to 
the data is possible, we should follow a similar 
program for the color-singlet rate |[2|. This will 
allow us to resum subleading Sudakov logs. The 
existence of these logs was first pointed out in jllj , 
where it was observed that though the leading 
logs cancel in the color-singlet differential rate, 
they are present in the derivative of the rate. 
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